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We draw an analogy between the chiral extrapolation of lattice QCD calculations from large 
to small quark masses and the interpolation between the large mass (weak field) and small mass 
(strong field) limits of the Euler-Heisenberg QED effective action. In the latter case, where the 
exact answer is known, a simple extrapolation of a form analogous to those proposed for the QCD 
applications is shown to be surprisingly accurate over the entire parameter range. 



The challenge to find an accurate and reliable method of chiral extrapolation for hadronic properties calculated in 
lattice QCD at large quark mass is a matter of considerable current importance. While computer limitations mean 
that lattice simulations at physical quark masses are many years away, recent progress in chiral extrapolation suggests 
that it may well be possible to obtain accurate hadronic properties based on the calculations which will be possible 
with the next generation of supercomputers, available within just a few years, in the 10 Tera-flops range. 

Fundamental to this scheme is the development of extrapolation methods which incorporate the model independent 
constraints of chiral symmetry [jl]Q] , notably the leading non-analytic (LNA) behaviour of chiral perturbation theory 
HU , as well as the heavy quark limit § . Although these extrapolations are designed to match the leading behaviour 
in the extreme limits of small and large quark mass, there has been little guidance as to their reliability in the 
intermediate mass region. It is very unclear what precision to expect from such a simple extrapolation into the 
intermediate mass region, because the large mass expansion is presumably asymptotic, and the small mass limit 
has a log divergence plus finite corrections with a small radius of convergence. Here we attack this question from a 
novel direction by considering a remarkably close analogy between this problem and a well-known, exactly soluble 
system in effective field theory — the Euler-Heisenberg effective action . The Euler-Heisenberg system exhibits 
many of the features found in the QCD calculations: at small electron mass (equivalently, strong external field) there 
is a logarithmic branch point, while at large mass (equivalently, weak external field) one has an asymptotic series 
expansion in inverse powers of mass. In this Letter, we show that a simple two-parameter interpolation formula (of 
the form used in the context of chiral extrapolation) , which builds in the correct leading behaviour in both the small 
and large mass limits, yields an excellent approximation to the exact Euler-Heisenberg answer over the entire range 
of mass. We discuss possible consequences of this observation for the chiral extrapolation of lattice data. 

Effective field theory (EFT) plays an important role in modern theoretical physics Ij-jnJ. In pioneering work in 
the 1930's, Heisenberg and Euler ||, and Weisskopf 0, studied the quantum corrections to classical electrodynamics 
associated with vacuum polarization effects. Renormalization properties and a more formal "proper-time" version 
were later studied by Schwinger In modern language, they computed the low energy effective action for the 
electromagnetic field, to leading order in the derivative expansion, by integrating out the electron degrees of freedom 
in the presence of a constant background electromagnetic field. This one-loop effective action can be expressed as 



where Ip — Y (d v + ieA v ), and A v is the fixed classical gauge potential with field strength tensor = d^A„ — dyA^. 
As shown in |3]-|8) , this effective action can be computed in a simple closed form when the background field strength 
F^ u is constant. For simplicity, we consider the case when the background is a constant magnetic field of strength B 
(and we choose eB to be positive). Then the exact, renormalizcd, onc-fcrmion-loop effective action has the following 
integral representation: 
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The j term is a subtraction of the zero field (B = 0) effective action, while the | subtraction corresponds to a 
logarithmically divergent charge renormalization M . 

We stress that Eq. (^) is an exact, non-perturbative result. However, it can of course be expanded in two obvious 
limits. In the large mass limit, m 2 3> eB (which is equivalently the weak field limit), it is straightforward to develop 
an (asymptotic) expansion of this integral: 
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Here the Bm are the Bernoulli numbers The large mass expansion, Eq. (||), of the effective action has the 

standard form, S = m 4 J2 n 01 a l° w energy effective action ||[l0]], where the higher-dimension operators O'"' (of 
dimension n) are balanced by n powers of the mass scale m, below which the low energy effective action is meaningful. 
In this case, "low energy" means that the cyclotron energy scale ^ is much smaller than the energy scale set by 
the electron mass m. That is, ^| <C 1. An alternate perspective on the large mass expansion is as a perturbative 
expansion in powers of the coupling e, with the n th power of e being associated with a one-fermion-loop diagram with 
n external photon lines (the divergent 0(e 2 B 2 ) self-energy term is not included, as it contributes to the bare action 
by charge renormalization Q). We note that, as a consequence of charge conjugation (Furry's theorem), only even 
powers of appear in the perturbative expansion of Eq. (|^) . It is interesting to note that the series expansion of 

Eq. @) is divergent, because the Bernoulli numbers grow factorially as Bin ~ 2(— 1)" +1 for large n, consistent 

with very general results for perturbation theory [ fbl|]l5| ]. It is in fact an asymptotic series, and the proper-time 
integral representation in Eq. (^) is just the straightforward Borel sum |l(| of this asymptotic series [l7| ]. 

The large mass limit may equivalently be characterized by the relevant length scales: the electron Compton wave- 
length A e = — , and the cyclotron radius ("magnetic length") Xb = -J=g- in terms of these length scales, the large 
mass limit corresponds to the situation where the electron Compton wavelength is much smaller than the cyclotron 
radius: A e < Ab. 

Since the Euler-Heisenberg system is exactly soluble, we can also use the exact integral representation (||) to study 
the small mass, or strong field, limit where m 2 <C eB. In terms of the length scales, in this limit the electron 
Compton wavelength is much greater than the cyclotron radius: A e 3> A^. Then, from Eq. (^), one finds (using 
results in Ref. fill): 
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Note that the coefficient, — fjfr, of the leading term, the log ^ term, is fixed by the (one-loop) QED beta function 
In (|), 7 is Euler's constant, and ((s) is the Riemann zeta function @. Note that C'(-l) ~ -0.165421. 
It is instructive to contrast this small mass expansion, Eq. (Q), with the large mass expansion, Eq. (j^). In the small 
mass limit, analogous to the chiral limit in QCD, we see the appearance of logarithmic terms, analogous to the "chiral 
logs" of QCD. In addition, note that both even and odd powers of ^ appear in the small mass expansion, Eq. (Q). 
On the other hand, in the large mass expansion, Eq. (^), there are no non-analytic log terms, and only even powers 
of ^| appear. So, we see that the one-loop Euler-Heisenberg effective action, which is given by the exact integral 
representation (|^), has two very different expansions in the two limits of large and small electron mass. The transition 
between these two extreme regions is governed by whether the electron Compton wavelength, A e , is larger or smaller 
than the cyclotron radius, A^. In Fig. [I] we plot the exact Euler-Heisenberg effective action, Eq. (||), with an overall 
factor of — e 8 ^| removed, as a function of and compare it to the leading large mass term — ^ (^) from Eq. (|^), 
and to the leading small mass terms | log ^ -f 0.763969 from Eq. (Q). From this figure it is clear that these leading 
terms accurately capture the extreme behaviours of the exact result, but do not interpolate in the intermediate region 
where the scales are comparable. 
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FIG. 1. Comparison between the exact action (solid curve) for the Euler-Heisenberg model and the leading terms in the 
expansions about the weak (dashed curve) and strong field (dash-dot curve) limits. Note that m 2 is measured in units of eB. 

Having reviewed these pertinent aspects of the Euler-Heisenberg effective action, we now turn to what appears at 
first glance to be a completely different problem: the calculation of hadron properties as a function of quark mass, or 
through the Gell-Mann-Oakes-Renner relation (m 2 oc m q ), pion mass. Chiral perturbation theory permits a rigorous 
expansion of hadron properties about the chiral limit, where — > 0. For example, for the nucleon charge radius one 
finds M 
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where ± refers to the proton or neutron respectively. (Here fj, just sets the scale against which the pion mass is 
measured. It is arbitrary in the sense that a change in fj, is equivalent to a change in the constant term, ci.) Note 
that the charge radius diverges logarithmically in the chiral limit, with a model independent coefficient 
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On the other hand, in the large m w limit, heavy quark effective theory suggests that the charge radius should decrease 
as 
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plus higher inverse powers of m 2 . (In the heavy quark limit one has essentially a Coulombic problem and the charge 
radius is proportional to the Bohr radius which goes as l/m q and hence l/m w .) 

As discussed at length in Ref. ||, current lattice data for charge radii are confined to pion masses greater than 
600 MeV. The corresponding pion Compton wavelength, A^ , is then smaller than the calculated charge radius, which 
we may take as an indication of the size, R, of the source of the pion field. The lattice data shows only a very slow 
variation of {r 2 )E in the mass range where the lattice calculations have been made, with no indication of a chiral log. 
Yet, in order to compare with the physical charge radii one must extrapolate these lattice results to the chiral regime 
where X v ^> R and the chiral log is important. This is the challenge of chiral extrapolation. 

We wish to draw an analogy between the Eulcr Heiscnbcrg system discussed above and this system. In this analogy, 
the pion Compton wavelength, A^, plays the role of the electron Compton wavelength, A e , and the source size, i?, 
plays the role of the magnetic cyclotron radius, Xb, (equivalently, the mass scale /i 2 plays the role of the magnetic field 
strength eB). The chiral perturbation theory expansion of Eq. (||), where A^ ^> R, is analogous to the leading terms 
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in the small mass expansion of Eq. (|J), where A e 3> Xb- The heavy quark effective theory result presented in Eq. (|?|), 
where <C R 7 is similarly analogous to the leading term in the large mass expansion in Eq. (|^) where A e C A^. 

In the QCD context, following earlier studies of magnetic moments where it was found that a simple Pade 
approximant was able to describe the mass dependence arising in a particular chiral quark model, Hackett- Jones 
et al. J3| extrapolated the lattice data from m\ > 0.4 GeV 2 to m\ — 0.02 GeV 2 (the physical point) using an 
interpolating formula which was chosen as the simplest two-parameter form consistent with the constraints imposed 
by the extreme behaviours in the large and small pion mass limits, Eq. (Q) and Eq. (||) respectively. (Recall that xn 
is model independent, and note that the data could constrain no more than two parameters.) In the light of later 
experience |IJ , we choose to use a slightly modified argument in the chiral log: 
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Here, rather than being arbitrary, fi assumes physical significance as the scale above which the chiral log is suppressed 
— of course, Eq. (|J) preserves the correct behaviour in the chiral limit. From experience with moments of structure 
functions, magnetic moments and hadron masses, this scale is expected to be /i ~ 500 MeV. As the lattice data is 
not yet able to constrain /i, we simply fix it to 500 MeV and adjust only c\ and C2- Figure || shows the resulting fit 
to the proton charge radius and the corresponding extrapolation to the physical pion mass. As discussed in this 
chiral extrapolation fit is closer to the physical value than a naive linear fit through the lattice data. However, in the 
absence of lattice data at lower quark masses, it is difficult to be more precise about the quality of the fit. 
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FIG. 2. Fit to the lattice QCD data for the square of the proton charge radius as a function of pion mass squared, using 
Eq. (^) . The extrapolated value at the physical pion mass (indicated by the vertical dotted line) is shown by the solid dot with 
the large error bar, while the star indicates the experimentally observed value. 

In view of the close parallel between this hadronic problem and the Euler-Heisenberg system in QED, we return to 
the Euler-Heisenberg system, where we can be much more quantitative concerning the accuracy of an interpolating 
fit. We ask the following question. Suppose that we did not know the exact integral representation answer d) for the 
effective action, but that we did know the leading terms in each of the extreme large and small mass limits. Would it 
then be possible to find a simple two-parameter interpolating formula, analogous to (||), that connected the extreme 
limits in a smooth manner? And if so, how accurate would such an interpolating formula be in the intermediate 
region? 

The leading terms are determined as follows. In the large mass limit, this is the first term, (ffp) , in @i 

corresponding to the first nonlinear correction to classical electrodynamics, whose coefficient comes from the one- 
fermion loop with four external photon lines, a straightforward perturbative calculation. In the small mass limit, 
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the leading term in (Q) is the logarithmic term, — logjg-, whose coefficient is fixed by the one- loop QED 

beta function [Q. Motivated by the interpolation formula, Eq. (||), which was used in the QCD case, we propose the 
following interpolating function for the effective action 



'interpolating 



8tt 2 



dt + | log 



m 2 +eB 



(9) 



This interpolating formula has the correct leading behaviour in both the large and small m limits. Figure ^ shows a 
comparison of the fit obtained with this form by adjusting the two parameters d± and di (dash-dot curve) with the 
exact result (solid curve). Our best fit was obtained with parameter values: d\ — 0.7059, and c?2 = 1.5541. Figure 
3 also shows the percentage difference between the exact result and approximate expressions (dashed line). (Note 
that m 2 is expressed in units of eB.) Over the entire range of ^g-, the interpolating function is within 10% of the 
exact answer. Such precision is very surprising when we recall that the Euler-Heisenberg effective action has the 
problems (shared by the analogous QCD calculations) that the large mass expansion is asymptotic and the small 
mass expansion has a log divergence and a small radius of convergence. 
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FIG. 3. Comparison between the exact expression for the action in the Euler-Heisenberg model (solid line) and the interpo- 
lating approximation given in Eq. ^ which builds in the correct chiral and heavy quark limits (dot-dashed line). Note that the 
agreement is so good that it is difficult to distinguish between the two curves on this scale. The percentage difference between 
the two is indicated by the dashed line. 
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In summary, the Euler-Heisenberg system presents a problem which exhibits many of the mathematical compli- 
cations of the chiral extrapolation problem in QCD, yet it is exactly soluble. By carefully respecting both the high 
and low mass limits of the exact solution, we showed how to construct a simple formula which reproduced the exact 
solution over the entire parameter range with surprisingly good accuracy. Of course, in the Euler-Heisenberg case we 
have the advantage of fitting the exact function over the entire mass range, while in QCD we have to extrapolate from 
large quark mass (where lattice data is available) to the chiral limit. Nevertheless, the fact that the the mathematical 
structure of the two problems is identical, combined with the success achieved in the Eulcr-Hcisenberg problem, gives 
us considerable confidence that a similar level of accuracy may be obtainable for QCD. It is therefore extremely 
encouraging that the chiral extrapolation of even the present crude lattice data at very large quark masses yields a 
physical proton charge radius within one standard deviation of the experimental value. Even more important, this 
result lends enormous impetus to the quest for new lattice data at lower quark mass which will better constrain the 
chiral extrapolation. It suggests that the next generation of supercomputers (available within 2-3 years) may well 
provide sufficient information that, in combination with these chiral extrapolation techniques, one should be able to 
calculate accurate hadron properties at the physical quark mass. 
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